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We describe the application of the continuous wavelet transform to calculation of the Green 
functions in quantum field theory: scalar (j)'^ theory, quantum electrodynamics, quantum chro- 
modynamics. The method of continuous wavelet transform in quantum field theory presented 
in [1] for the scalar 0* theory, consists in substitution of the local fields (/>(a;) by those depen- 
dent on both the position x and the resolution a. The substitution of the action S[4){x)] by 
the action S[(l)a{x)] makes the local theory into nonlocal one, and implies the causality condi- 
tions related to the scale a, the region causality [2]. These conditions make the Green functions 
G{xi, ai, . . . ,Xn,an) = {(l>ai{xi) . . . (f)a,^{xn)) finite for any given set of regions by means of an effec- 
tive cutoff scale A — min(ai, . . . , a„). 

PACS numbers: 03.70.-(-k, ll.lO.-z 
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I. INTRODUCTION 

The fundamental problem of quantum field theory and 
statistical mechanics is the problem of divergences of 
Feynman integrals emerging in Green functions. The 
formal infinities appearing in perturbation expansion of 
Feynman integrals are tackled with different regulariza- 
tion methods, from Pauli-Villars regularization to renor- 
malization methods for gauge theories, see e.g. [3] for a 
review. A special class of regularizations are the lattice 
regularizations tailored for the precise numerical simula- 
tions in gauge theories [4, 5]. 

There are a few basic ideas connected with those reg- 
ularizations. First, certain minimal scale L = ^, where 
A is the cutoff momentum, is introduced into the the- 
ory, with all the fields (/)(x) being substituted by their 
Fourier transforms truncated at momentum A. The phys- 
ical quantities are then demanded to be independent on 
the rescaling of the cut-off parameter A. The second 
thing is the Kadanoff blocking procedure [6], which av- 
erages the small-scale fluctuations up to a certain scale - 
this makes a kind of effective interaction. 

Physically all these methods imply to the self- 
similarity assumption: blocks interact to each other sim- 
ilarly to the sub-blocks [7] . Similarly, but not necessarily 
having the same interaction strength - the latter can be 
dependent on scale A = A (a). However there is no place 
for such dependence if the fields are described solely in 
terms of their Fourier transform - except for the cutoff 
momentum dependence. The latter representation, being 
based on the representation of the translation group, is 
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rather restrictive: it determines the effective interaction 
of all fluctuations up to a certain scale, but says nothing 
about the interaction of the fluctuations at a given scale. 

That is why, the functional methods capable of taking 
into account the interaction at specific scale are required. 
Wavelet analysis, being the multiscale alternative to the 
Fourier transform, emerged in geophysics [8] , is the most 
known of such methods. Its application to quantum field 
theory have been suggested by many authors [9-13] The 
other side of the problem is that the quantum nature 
of the fields considered in quantum field theory is con- 
strained by the Heisenberg uncertainty principle. To lo- 
calize a particle in an interval Ax the measuring device 
requests a momentum transfer of order Ap ~ h/Ax. If 
Ax is too small the field (t){x) at a fixed point x has no 
experimentally verifiable meaning. What is meaningful, 
is the vacuum expectation of product of fields in certain 
region centered around x, the width of which (Ax) is 
constrained by the experimental conditions of the mea- 
surement [1]. That is why, at least from physical point 
of view, any such field should be designated by the reso- 
lution of observation (jjAxix). 

In present paper we exploit the observation that quan- 
tum field theory models, which yield divergent Feyn- 
man graphs, can be studied analytically if we project 
original fields (j){x) into the fields (j)a{x), subscribed 
by the scale of measurement a. The Green functions 
{(j)ai{xi) . . . 4>an{xn)) bccome finite under certain causal- 
ity assumptions, which stand for the fact that any n-point 
correlation function can be dependent only on space-time 
regions, rather than points, and thus cannot be infinite 
[1]. These Green function describe the effect of prop- 
agation of a perturbation from a region of size a, cen- 
tered at a point x, to a region of size a', centered at 
a point x' . The standard quantum field theory models 
can be reformulated by expressing the (point-dependent) 
local fields (j>{x), the distributions, in terms of the region- 
dependent fields (j)a{x). The integration over all scales a 



will of course drive us back to the known divergent re- 
sults, but the physical observables are always those mea- 
sured with finite resolution, and their correlations are 
always finite. Therefore the idea of wavelet transform of 
quantum fields, which will be considered below, is very 
similar to the idea of renormalization group (This simi- 
larity, being studied in the lattice framework [12, 14], is 
out of coverage of the present paper). 

The remainder of this paper is organized as follows. In 
Section II we recall the basics of the continuous wavelet 
transform and its application to the multiresolution anal- 
ysis of quantum fields. The definitions of scale-dependent 
fields and Green functions, the modifications of the Feyn- 
man diagram technique are presented. The 0* scalar 
field model examples of calculations are given. Section 
III considers the case of operator- valued scale-dependent 
fields. The operator ordering and commutation relations 
are presented. The relations between the theory of scale- 
dependent fields in Euclidean and Minkowski spaces are 
discussed. Section IV presents the examples of calcu- 
lation of one- loop Feynman graphs in QED and QCD. 
The Conclusion gives a few remarks on the perspectives 
and applicability of the multiscale field theory approach 
based on continuous wavelet transform. 



II. CONTINUOUS WAVELET TRANSFORM IN 
QUANTUM FIELD THEORY 

A. Basics of the continuous wavelet transform 

Let "H be a Hilbert space of states for a quantum field 
\(j)). Let G be a locally compact Lie group acting tran- 
sitively on H, with dfi(h'),h' G G being a left-invariant 
measure on G. Then, similarly to representation of a 
vector |(/)) in a Hilbert space of states "H as a linear 
combination of an eigenvectors of momentum operator 
\(f)) = J \p)dp{p\(j)), any |(/)) € % can be decomposed with 
respect to a representation U{v) of G in "H [15, 16]: 

\'t>) = ^ f U{iy)\g)d^i{:.){g\U*{,^M), (1) 

^g JG 

where \g) G "H is referred to as an admissible vector, or 
basic wavelet, satisfying the admissibility condition 



G ^ 

'"h¥JG 



\{g\U{v)\g)\^d^l{v)<^. 



The coefficients {g\U*{v)\(t>) are referred to as wavelet 
coefficients. 

If the group G is abelian, the wavelet transform (1) 
with G : x' = X -\- h' coincides with Fourier transform. 

The next to the abelian group is the group of the affine 
transformations of the Euclidean space W^ 

G:x'^ aR{e)x + fe, .x, 6 G M'^, a G M+, G SO{d), (2) 

where R{0) is the rotation matrix. We define unitary 
representation of the affine transform (2) with respect to 



the basic wavelet g{x) as follows: 

C/(a,6,%(x) = lgfi?-i(0)^— ^). (3) 

(We use L^ norm [17, 18] instead of usual IP' to keep the 
physical dimension of wavelet coefficients equal to the 
dimension of the original fields) . 

Thus the wavelet coefficients of the function 0(x) G 
L'^iW^) with respect to the basic wavelet g{x) in Eu- 
clidean space R'' can be written as 



<Pa.m^ I -^9{R-'{0)'^^\{x)d''x. (4) 

The wavelet coefficients (4) represent the result of the 
measurement of function 4'{x) at the point b at the scale 
a with an aperture function g rotated by the anglc(s) 6 
[19]. 

The function 4){x) can be reconstructed from its 
wavelet coefficients (4) using the formula (1): 

^(^) -^ l-.a {R-\ef—^) <^,,(fe)^^dM^) 

Cg J a'^ \ a J a 

(5) 
The normalization constant Gg is readily evaluated using 
Fourier transform: 



Go 



,da 



\g{aR-\9)k)\^-d^l{e) = / |.9(fc)| 



.d'^k 



< oo. 



For isotropic wavelets 
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\~giak)\ 



,da 



a 

d/2 
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d'^k 



where Sd = YUh) ^^ ^^^ area, of unit sphere in R''. 



B. Resolution-dependent fields 

If the ordinary quantum field theory defines the field 
function (t>{x) as a scalar product of the state vector of 
the system and the state vector corresponding to the lo- 
calization at the point x: 

<i>{x) = (x|0), (6) 

the modified theory [1, 20] should respect the resolu- 
tion of the measuring equipment. Namely, we define the 
resolution- dependent fields 



(paeix) = {x,0,a;i 



(7) 



also referred to as the scale components of (j), where 
{x,9,a;g\ is the bra-vector corresponding to localization 
of the measuring device around the point x with the spa- 
tial resolution a and the orientation 9 G SO{d); g labels 
the apparatus function of the equipment, an aperture [19]. 



The field theory of extended objects with the basis g de- 
fined on the spin variables was considered in [21, 22]. 

The introduction of resolution into the definition of the 
field function has a clear physical interpretation. If the 
particle, described by the field 4>{x), have been initially 
prepared in the interval (x— ^, a;+ ■^), the probability 
of registering this particle in this interval is generally less 
than unity: for the probability of registration depends on 
the strength of interaction and the ratio of typical scales 
of the measured particle and the measuring equipment. 
The maximum probability of registering an object of typ- 
ical scale Ax by the equipment with typical resolution a 
is achieved when these two parameters are comparable. 
For this reason the probability of registering an electron 
by visual range photon scattering is much higher than 
by that of long radio-frequency waves. As mathematical 
generalization, we should say that if a measuring equip- 
ment with a given spatial resolution a fails to register 
an object, prepared on spatial interval of width Ax with 
certainty, then tuning the equipment to all possible res- 
olutions a' would lead to the registration. This certifies 
the fact of the existence of the measured object. 

In terms of the resolution-dependent field (7) the unit 
probability of registering the object (p anywhere in space 
at any resolution and any orientation of the measuring 
device is expressed by normalization 



\(i>a.0{x)\'^dn{a,0,x) = 1, 



(8) 



where dfj,{a, 9, x) is an invariant measure on R+ x SO{d) x 
R'^, which depends on the position x, the resolution a, 
and the orientation 9 of the aperture g. 

If the measuring equipment has the resolution A, i.e. 
all states {g;a > A,x\(j)) are registered with significant 
probability, but those with a < A are not, the regu- 
larization of the field theory in momentum space, with 
the cutoff momentum A ~ 2tt/A corresponds to the UV- 
regularized functions 



0(^)(x)= 



1 

a, 



>A 



{x\g;a,b)dfi{a,b){g;a,b\(p). (9) 



The regularized n-point Green functions are 
g(^)(xi, . . . ,x„) = (0(^)(xi), . . . ,0(^)(x„)),. 

However, the momentum cutoff is merely a technical 
trick: the physical analysis, performed by renormaliza- 
tion group method [3, 23, 24], demands the independence 
of physical results from the cutoff at A — > oo. 



C. Scalar field example 

To illustrate the method, following [1, 20], we start 
with Euclidean scalar field theory. The widely known 
example which fairly illustrates the problem is the 0^ 
interaction model in M"^, see e.g. [3, 25], determined by 
the generating functional 

Vt.[j] ^^ /•,-/'i^-[^W)^+^*^+A<A^-Hp^, (10) 

where A/" is a formal normalization constant. The con- 
nected Green functions are given by variational deriva- 
tives of the generating functional: 



Q{n) ^ 



(5"lnTy[J] 



SJ{xi) . . .6J{Xn) 



(11) 



J=0 



In statistical sense these functions have the meaning of 
the n-point correlation functions [26]. The divergences 
of Fcynman graphs in the perturbation expansion of the 
Green functions (11) with respect to the coupling con- 
stant A emerge at coinciding arguments Xi — x^. For 
instance, the bare two-point correlation function 



G^^\x-y) = 



d'^p e 



-ip(x-y) 



(277)^^ p2 _^_ ^2 



(12) 



is divergent at x — y for d> 2. 

For simplicity let us assume the basic wavelet g to 
be isotropic, i.e. we can drop the rotation matrix R{9). 
Substitution of the continuous wavelet transform (5) into 
field theory (10) gives the generating functional for the 
scale-dependent fields 4>a{x) [20]: 



Ww[Ja] 



N I 2?(/)a (a;) exp -- / (pai{xi)D{ai,a2,xi 



X2)4'a2{x2) 



^ / K?;.::;i'!r0a,(xi)---^a.(a:4) 



daid'^xi da2d'^X2 da^d'^x^ da4d'^X4 



ai 



02 



03 



04 



daid'^xi da2d'^X2 



oi 



Ja{x)(j)a{x) 



0.2 

dad'^x 



(13) 



with D(ai, 02, xi — X2) and V^"^;.;;;£f denoting the wavelet 
images of the inverse propagator and that of the interac- 
tion potential. The Green functions for scale component 
fields are given by functional derivatives 

S^ In Ww[Ja] 



{(l>aiixi)-- ■(f)a„iXn))c = 



SJaAxi) ■■■SJa„{x„) 



,7=0 



Surely the integration in (13) over all scale variables 
/o°° ^ turns us back to the divergent theory (10). 

This is the point to restrict the functional integration 
in (13) only to the field configurations {4'aix)}a>A- The 
restriction is imposed at the level of the Feynman dia- 
gram technique. Indeed, applying the Fourier transform 



to the r.h.s. of (5,4) one yields 

I f^ da f df^k 



C, 



g JO 



-ikx y. 



(27r) 



{ak)(i>a{k), 



4>a{k) = g{ak)(j){k). 



Doing so, we have the foUowing modification of the Feyn- 
nian diagram technique [13]: 

• each field (t){k) will be substituted by the scale com- 
ponent (j){k) -^ (paik) = g{ak)(t>{k). 

• each integration in momentum variable is accom- 
panied by corresponding scale integration: 



d'^k 



d'^k da 



{2'kY C^ttY a 



• each interaction vertex is substituted by its wavelet 
transform; for the A^-th power interaction vertex 

N 

this gives multiplication by factor I I g(aiki). 



According to these rules the bare Green function in 
wavelet representation takes the form 



G'-^\ai,a2,p) = 



g{aip)g{-a2p) 



pA 



The finiteness of the loop integrals is provided by the 
following rule: there should be no scales ai in internal 
lines smaller than the minimal scale of all external lines. 
Therefore the integration in a^ variables is performed 
from the minimal scale of all external lines up to the 
infinity. 

To understand how the method works one can look at 
the one-loop contributions to the two-point Green func- 
tion G^^^ai,a2,p) shown in Fig. la., and to the vertex 
shown in Fig. lb. The best choice of the wavelet function 



a) 



'¥■ 




^ 



b)' 






X ^X^)^ 



^8 3 



of [20], is written as 



1 



rf(Am)= — 



q Ja,^a,>A {27rr 



d'^q l.9(«3g)|^|g(-Q4g)|^ ^03 duj 

as a4 
x'^^^dx 



bdm ' 



f^{Amx 



(27r)'i Jo •' '^^""■"' a;2-Kl ' 

where the integration over the scale variables resulted in 
the effective cutoff function 



fix)^ 



9 Jx 



\~a{a)f-. f{0)^l. (14) 

a 



In the one-loop contribution to the vertex, shown in 
Fig. lb, the value of the loop integral is 



X, 



A2 



1 



A .p{qA)fHiq-s)A) 

2 (27r)'' J (277)^^ [g2 + m2] [{q - s)^ + m 



(15) 



where s = pi+p2, A = Tam{ai, 02, 03,04). The integral 
(15) can be calculated by symmetrization of loop mo- 
menta <?—>'/+§ in Fig. lb, introducing dimensionless 
variable y = q/s, after a simple algebra we get 



\2 c „d— 4 f-TT coo 



(27r 



\2d 



dyy 



d-3. 



P (^As^y^ + ycos0 + l^ P (As^y^ - ycosO + \ 



v'+i+- 



COS( 



— cost 



where 9 is the angle between the loop momentum q and 
the total momentum s. For the simple choice of the basic 
wavelet gi [1, 20] 



9{x) 



o-x^/2 



g{k) = ike^^^'^ 



(27r)''/2 
in four dimensions we get a finite result 

-4a4e2"'Eii(2a2)_^2a2 ^ 



lim X4^{a ) = 



-2a 



s2^4m 



647r2a4 
2a2e2"'Eii(2a2)], 



depending on is dimensionless scale factor a = Am, where 
A is the minimal scale of all external lines. 

These results display an evident fact that for the mas- 
sive scalar field all length scales are to be measured in 
the units of inverse mass. 



+ permutations + ... 

FIG. 1: Feynman diagrams for the Green functions G^^' and 
G' ' for the resolution-dependent fields. Redrawn from [1] 

g{x) would be the apparatus function of the measuring 
device, however any well localized function with g{0) ~ 
will suit. The tadpole integral, to keep with the notation 



III. CAUSALITY AND COMMUTATION 
RELATIONS 

A. Operator ordering 

Up to now we have considered the calculation of the 
Feynman diagrams for the scale-dependent fields (f>a.- (x) 



treated as c-valued functions. In quantum field tlieory, 
adjusted to high, energy physics applications, the fields 
(j)a.-{x) are operator-valued functions. So, as it was al- 
ready emphasized in the context of the wavelet appli- 
cation to quantum chromodynamics [1, II], the operator 
ordering and the commutation relations are to be defined. 
In standard quantum field theory the operator ordering 
is performed according to the non-decreasing of the time 
argument in the product of the operator- valued functions 
acting on vacuum state 

A{ta)A{tn^i)...A{h)A{h)\0). 

t„>t„-l>...>t2>ti 

In the infinite momentum frame (IMF), which simpli- 
fies algebraic structure of the current algebra, the time- 
ordering is performed in the proper time argument x_|_ 
[27]. The quantization is performed by separating the 
Fourier transform of quantum fields into the positive- and 
the negative-frequency parts 



(b = ^+{x) + ^ (x). 



defined as follows 



d'^k 



^[e*'=^u+(fc)-^e-^'^^u-(fc)], (16) 



(2^) 

where the operators u^{k) — u{±k)9{ko) are subjected 
to canonical commutation relations 

[u+{k),u-{k')] = A{k,k'). 

In case of the scale-dependent fields, because of the pres- 
ence of the scale argument in new fields (j)a,riix), where 
a and t] label the size and the shape of the region cen- 
tered at X, the problem arises how to order the operators 
supported by different regions. This problem was solved 
in [1, 28] on the base of the region causality assumption 
[2]. If two regions {Ax,x) and {Ay,y) do not intersect 
the standard time ordering procedure is applied. Alter- 
natively, if one of the regions is inside another, see Fig. 2 
the operator standing for the bigger region acts on vac- 
uum first [I]. This causal ordering, drawn in Euclidean 
space, is presented in Fig. 2 below. The time ordering in 
Euclidean space, as an analytic continuation of time or- 
dering in Minkowski space, have been already considered 
in [29]. The diagram Fig. 2 shows space- like regions in 
Euclidean space. For Minkowski space corresponding di- 
agrams can be obtained by analytic continuation of the 
Euclidean ball of imaginary radius lA into Minkowski 
space, where we can restrict ourselves with forward light 
cone i > 0, \x\ < t. The disjoint events in Minkowski 
space are shown in Fig. 3. The correspondence to the 
other case of one Euclidean event inside another, shown 
in Fig. 2b, looks more complex after analytic continua- 
tion to Minkowski space. The forward light-cone part of 
such intersection is shown in Fig. 4. 

Phenomenologically, the principle 'the coarse acts first' 
is related to the definition of the measurement procedure. 



0- 



Y 




AY 




a) 



FIG. 2: Causal ordering of scale-dependent fields. Space-like 
regions are drawn in Euclidean space: a) The event regions 
do not intersect; b) Event X is inside the event Y. 




FIG. 3: Disjoint events in (f , a::) plane in Minkowski space 



possibly generalized, where the state of a part can be 
measured/affected only after and relatively to the state 
of the whole. The similar reason underlies the restriction 
on the scales in internal loops by the minimal scales of 
all external lines of the Feynman diagram: if we measure 
a quantum system from outside we cannot excite modes 
finer than the minimal available scale of measurement. 
Thus the functional integration over the trajectories in 




FIG. 4: Nontrivial intersection of two events X C Y in {t, x) 
plane in Minkowski space 

the space of square-integrable functions ^(^{x) is sub- 
stituted by functional integration over all causal paths, 
or tubes of all different thickness, in the space of scale- 
dependent functions V^ia {x) . Referring the reader to the 
original works of [2, 30] for the topological aspects of 
causal paths we ought mention that the Bogolioubov mi- 



crocausality condition holds for causal paths in the same 
way as it holds for trajectories [1]. It is also easy to 
show, that if the domain Y is inside the domain X the 
corresponding Green function is not singular at coincid- 
ing arguments - it is a projection from coarser scale to 
finer scale: 



G^^\ai,a2,bi-b2 = 0) 



d^p g{aip)g{-a2p) 
(27r)4 p^ + m^ ' 



-tp-O 



since \g{p)\ vanish at p — > oo. 



B. Commutation relations 

In case of wavelet transform the positive- and negative- 
frequency part operators (16) can be expressed using 
wavelet transform 

ufik) = 7^ / rfW -~g,iaM-\rj)k)ut^{k), 



from where we can set [1]: 



(17) 



[u+{k),u-ik')] (18) 



to ensure canonical commutation relations for 

[u+{k),u-{k')]. 



C. The Dyson-Schwinger equation 



Ordering in scale argument results in the modification 
of the Dyson-Schwinger equation in the theory of scale- 
dependent functions. Let G{x—y, Ux, ay) be the bare field 
propagator, describing the propagation of the field from 
the region (y, ay) to the region (x, a^)- Let Q{x—y, ax,ay) 
be the full propagator for the same regions. The Dyson- 
Schwinger equation relating the full propagator with the 
bare propagator is symbolically drawn in the diagram 



ay a^ ay a^ ay af^^a2 a^ 



(19) 



The integral equation depicted in diagram (19) can be 
written as 



g{x-y,ax,ay) ^ G{x-y,ax,ay)+ / / dxidx2G{x - X2,ax,a2)V{x2 - Xi,a2,ai)g{xi ''y,ai,ay), (20) 

J ai J a2 J 

I 

where the full vertex ^^ = 'P{x2 — xi, 02,01) denotes boson, or the self-energy diagram otherwise. The Fourier 
the vacuum polarization operator if G is the massless counterpart of the equation (20) can be written as 



Ga^,ay{p) = Ga^,ay{p) + / "^ 



Ga^.a2{p)'Pa2.ai{p)Qat,ay{p)- 

02 



r 



D. Wavelet transform in Minkowski space 



The straightforward application of wavelet transform 
(4), defined in Euclidean space R'', to the Minkowski 
space Mf would be the analytic continuation of the re- 
sults into the imaginary time 2:4 = ixq, making the Eu- 
clidean rotations into Lorentz boosts. The construction 
of such wavelets with respect to the representations of 
the Poincare group have been studied by several authors 
[31, 32]. From physical point of view there exists a sim- 
ple and an elegant way of making wavelet transform in 
Minkowski space. 



In quantum field theory problems related to relativis- 
tic particle collisions we can always change the coordinate 
frame to the co-moving frame of a relativistic projectile 
moving at utmost speed of light. Due to the Lorentz 
contraction of the projectile the longitudinal and the 
transversal degrees of freedom behave essentially differ- 
ent in such system. Without loss of generality we can 
always assume the projectile to move along the z axis. 

The Lorentz contraction, i.e., the hyperbolic rotation 
in (t, z) plane is determined by the hyperbolic rotation 
angle - the rapidity 77. The rotations in the transverse 
plane are not affected by the Lorentz contraction and are 



determined by the 5*0(2) rotation angle (/). If the problem 
is axially symmetric, the latter can be dropped. 

Therefore it is convenient to change from the space- 
time coordinates {t,x,y,z) to the light-cone coordinates 
{x+,x-,x,y): 



x± 



t±z 



x_L = {x,y). 



(21) 



This is the so-called infinite momentum frame. The IMF 
is not a Lorentzian system, but a limit of that at f — > c. 
The advantage of the coordinates (21) for the calcula- 
tions, say in QED or QCD, is significant simplification of 
the vacuum structure [27, 33]. The metrics in the light- 
cone coordinates becomes 



9f^iy 



^0 1 
10 
0-10 

VO -1/ 




The rotation matrix has a block-diagonal form 

/e'' 
p"** 
M{v,^)= 

Vo 

so that M-i(r/,</)) = M{-'q, -(/)). 

Wc can define the wavelet transform in light-cone co- 
ordinates in the same way as in Euclidean space using 
the representation of the affine group 

x' = aM{r],(j))x + b, 

defined in L^ norm as 



Uia,b,7j,^)u{x) ^ ^u { M-'iri, 



-b 



we have the definition of wavelet coefficients of a function 
f{x) with respect to the basic wavelet g as follows 



Wa,b,rj,4,[f] = 



2 1 

dxj^dx^d xj^— -X 
a* 



X 5 M- 1(77,0): 



/(x+,x_,x_l). 



In contrast to wavelet transform in Euclidean space, 
where the basic wavelet g can be defined globally on R'', 
the basic wavelet in Minkowski space is to be defined sep- 
arately in four domains impassible by Lorentz rotations: 



Ai : fc+ > 0, fc_ < 0;A2 : fc+ < 0, fc_ > 0; 
A3 : fc+ > 0, fc_ > 0;A4 : fc+ < 0, fc_ < 0, 



where k is wave vector, k± = "^"^ . Whence we have 
four separate wavelets in these four domains [34] : 



g,{x) = / e^'^-gik) 



d*k 



79^^4- (23) 

We assert the following definition of the Fourier trans- 
form in light cone coordinates: 



,f(x+,x^,x^) = J e'''--++'''+--'''^^\f{k^,k+,-k^)x 

dk^dk^d'^kj^ 



Substituting the Fourier images into the definition (22) 
we get 



W, 



ab-nd) I ^ 



,tk—b+-\-tk+b— — zk^bi 



/(fc_,fc+,ki) 



~g{ae^k-,ae-'^k+,aR-\q^)k^)^—^. (24) 



In Fourier space the relation between Fourier coefficients 
and wavelet coeflficients is therefore the same as in W^: 



Wa^^ik) = f{k)giaM-\r^,cl))k). 



^ I Similarly, the reconstruction formula is [35]: 



/(^)- 



u 

2 — 1 ^ 



27r 



dj] \ d(j) — 



E?^ r^^/^^<^r- / db+db^d'h^^9JM-\rj)LJ.)wi 
^ Cg, J -00 Jo Jo a JMf a* V a / 



A, 



(2^)4 



-Warj4,ik)g{ak^e^ ,ak+e ^,aR ((/))kj 



If the problem is axially syminetric the azimuthal part of 
integration (0) can be dropped. It is also convenient to 
use logarithmic scale a = e" to study different rapidity 
domains. 



E. Choice of the basic wavelet 

The choice of the basis of wavelet decomposition is 
always a subtle question, specially in quantum field the- 
ory. (The best choice, as it was already emphasized in [I], 
would be the apparatus function of a classical measuring 
device interacting with quantum system.) Some basis 
is always tacitly assumed. Even describing the mass- 
less photons, which are not localized anywhere, by plane 
waves, the possibility of photon registration by photo- 
multiplier implies its interaction with electron, and hence 
some scale and some localization. 

If the continuous wavelet transform is used in place of 
the Fourier transform the choice of the basic function is 
constrained by the feasibility of the analytical integration 
in Feynman diagrams on one hand, and by the possibility 
to understand this basic function as a localized (quasi) 
particle. The latter being claimed by some authors to be 
important for Minkowski space [36], seems questionable 
for Euclidean space calculations. If the wavelet trans- 
form is performed on a lattice there is a bias that only 
the similarity properties are important, rather than the 
shape of wavelet [12, 14]. The question whether or not 
the basic wavelet should satisfy some equation of motion 
is still not clear. 

To justify our choice of the derivatives of the Gaussian 
as the basic wavelets, we present the following heuristic 
consideration, inferred from the coherent states theory 
[32]. Let us introduce a localized wave packet in Fourier 
space 



~g{t,k) 



_ -itk-k^/2 



(25) 



If the wave packet is considered in Minkowski space, then 
k^ = can be assumed for the photon and the whole 
equation (25) turns to be a plane wave. Otherwise it is a 
localized wave. If t is time the packet (25) is a gaussian 
wave packet at initial time t — 0. At finite but small 
instants of time the wavepacket can be approximated by 
its Taylor expansion 



git, k) = Uk) + yM^) + 2Mk) + 0{f'), 
where the expansion coefficients 






are responsible for the shape of the packet at the time 
scales at which 1, 2 or ri excitations are significant. 
Clearly gn{x) are the excitations of the harmonic oscilla- 
tor, with gi being the first excitation, see Fig. 5. 




FIG. 5: Graph of gi wavelet: gi{x) — —xe ^ 1^ 



IV. GAUGE THEORIES 

A. QED 

Quantum electrodynamics is the simplest case of gauge 
theory. The local C/(l) invariance of the fcrmion field 



V'(x) 



-xeK^x) 



i/'(x) 



is accompanied by the gradient invariance of the vector 
field Ap(x), the electromagnetic field. 



A^,{x) ^ A^{x) ^ ^^,!^{x), 



(26) 



to keep the total action S'(V',V':^) — J Ld*x invariant 
under the local U{1) transform generated by A(x). The 
interaction of the charged fermion field -0 with electro- 
magnetic field Afi is introduced by substitution of ordi- 
nary derivatives d^ to covariant derivatives 



Da = d„ 



°^\{x), 



with e being the charge of the fcrmion. 

The Lagrangian of QED has the (Euclidean) form 



1 



1 



L = i,{x){lp + im)ij{x) + -F^.^^''+ — (a^A^)' , 



2a 



gauge fixing 
with Ff,^ = df,A^ - d^A^ (27) 

being the field strength tensor of the electromagnetic field 
A, and "ffj, being the Dirac 7-matrices. 

The straightforward application of the Feynman inte- 
gral to the gauge theory with the Lagrangian (27) would 
be inefficient for the integration over the field A{x) would 
contain an infinite set of physically equivalent field con- 
figurations. For this purpose the gauge fixing, which re- 
stricts the integration only to gauge-nonequivalent field 
configurations was introduced by Faddeev and Popov 
[37]. 



Quantum electrodynamics is the most firmly estab- 
lished and most verified field theory model in physics of 
elementary particles. The probability amplitude of scat- 
tering obtained at a tree-level are in fairly good agree- 
ment with classical scattering theory. Starting from one 
loop level the Feynman integrals are formally divergent, 
and the physical results are derived using the renormal- 
ization invariance of QED. The most accurate tests for 
the renormalized calculations of the electron-photon in- 
teraction are the Lamb shift of the Hydrogen-like ion 
energy levels and the anomalous magnetic momentum of 
the electron [38-41]. 

In one-loop approximation the radiation corrections 
in QED come from three primitive Feynman graphs: 
fcrmion self-energy S(p), vacuum polarization operator 
n^^(p), and the vertex function Tp{p,q). In Euclidean 
space the equations for the above three graphs have the 
form: 



Electron self-energy 



S(p) 






(28) 



gives the corrections to the bare electron mass ttlq related 
to irradiation of virtual photons. 
Vacuum polarization diagram 



""-^^^^"''/(^'PtT.— ^ 



1 



-lu- 



(29) 



could be expected to give the nonzero corrections to the 
photon mass, but due to gauge invariance the photon 
mass remains zero, instead the one loop contribution 
(29) renormalizes the electron charge at small distances, 
therefore modifies the Coulomb interaction by screening 
the bare electron charge cq by virtual electron-positron 
pairs polarizing the vacuum at small distances. This di- 
agram contributes to the Lamb shift of the atom energy 
levels. 

One-loop vertex function 



T^pip,q) = -«e' 



rfV 



1 



1 



7p- 



(2^)4"^ + / + m"^/ + ^ + 



(30) 

determines the anomalous magnetic moment of the elec- 
tron. 

All three integrals (28,29,30) are divergent. Their 
evaluation involves regularization procedures. The most 
common is the dimensional regularization with all inte- 
grals taken in formal d — 2uj dimension with physical 
value u! ^ 2. In this way the divergences come as poles 
in e = 2 - w, see e.g. [23, 25, 42]. 



B. One loop corrections in wavelet-based theory 



similar to that of scalar theory (15). The evaluation of 
the one-loop radiative corrections for the scale-dependent 
fields give finite results by construction with no regular- 
ization procedure required. 

a. Electron self-energy. For the scale components of 
the electron self-energy diagram, we get 



S(-^)(p) 
9iap)9{-a'p) 



d^q FA{p,q)lt, 



It^ 



(2^)4 



(31) 
where A is the minimal scale of two external lines shown 
in Fig. 6: ^ = min(a, a'). The regularizing function 
Fa{p, q) is the result of integration over the scales of two 
internal lines. For the isotropic basic wavelet g the regu- 



p,a 



p/2+q 



X 



p/2-q 



> 



p,a 



FIG. 6: Electron self-energy diagram 

larizing function is given by (14): 

Fa{p, q) = .f{A{p/2 - q))f{A{p/2 + q)). (32) 

Introducing the dimensionless variable y — q/|p|, after 
straightforward algebra, we can perform the integration 
in Euclidean space: 



E(^)(p) 



d'y 



Kap)g{-a-'p) J (2^)4 

p + Am- 2|p|^ 



FA{pAp\y)^ (33) 



y^ + i 



ycostf- ^ 



r 



y cos 9] 



where 9 is the Euclidean angle between the p and the q 
directions. In high energy limit, p^ 3> 4?n^, the contribu- 
tion of last term in the numerator of (33) vanishes for the 
symmetry, and the equation (33) can be easily integrated 
in angle variable (CI): 



S(-^)(p) 
5(ap)5(-a'rt 

Ri{p) 



47r2 



Ri (p) (p + 4to) where: 



dyyFAip, \p\y) 



l3Hy) 



m = y+^- 



The integral i?i (p) is finite for any wavelet cutoff function 
(14). For the gi wavelet we get 



The evaluation of Feynman diagrams with fermions 
and gauge fields in wavelet-based Euclidean theory is 



Riip) = 



- ,-A'p' 



dyye 



-iA-'p-'y^ 



1- Wl 



1 
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After a simple algebra this gives 
^I(P) = ^:4^(2^VEii(AV) - 4AVEii(2AV) 

(34) 

6. Vacuum, polarization diagram The vacuum polar- 
ization diagram in quantum electrodynamics of scale- 
dependent fields is obtained by integration over the scale 
variables of the fermion loop shown in Fig. 7: 



2 / d'^'l 



'l{ap)g[-a'p) J (27r)4 

Sp(7m(^ +^/2 ^ rn)j,{^ - fl2 - m)) 
\{q + p/2)2 + m^] \{q - pjlf + m^] 



[(g+|)2+m2] [(g- 1)2+^2] 

(35) 

Similarly to previous diagram we use the dimensionless 

q+p/2 



p,a 




q-p/2 



p,a 



FIG. 7: Vacuum polarization diagram in (Euclidean) scale- 
dependent QED 

variable y and integrate over the angle variable. The mo- 
mentum integration in equation (35) is straightforward: 
having expressed all momenta in units of electron mass 
TO, we express the loop momentum in terms of the photon 
momentum and perform the integration over the angle 
variable: 



where A and B depend only on the modulus of y, but 
not on the direction, and can be expressed in terms of 
angle integrals (CI). 

Finally, writing the polarization operator as a sum of 
transversal and longitudinal parts, we have the equations 



l{ap)g{-a'p) 



A^) 



where 






Pt^Pv \ [A) 
„2 ) ^T 



X 



p. 

(A) PlJ-PiJ 



p^ 



2 2 



dyyFA{mp,mpy) 



(36) 



1- 

V ^ 

5 



16 ' w 



1 
2^ 






V 



p2 pi 



2 2/2 



7r[ = -—^m?p'^ j dyyFA{mp,mpy) 



2y^ 



-4y 



1 - 



16 ' » 



1 

2^ 



2y^ 



2yMi 



M) 



M^ 



g2^2p2 ,00 



dyyFA{mp,mpy) 



L^yi 



16 



1 
2^ 



■2y' 



\ \ ( 



V 



p^ p^ 8 



n 



(A) 



~, .~i ,v = ^(»TtV) / dyyFA{mp,mpy) X 

g{ap)g[-a'p) tt^ J„ 



X I d9 sin 
/o 



2y,y.^^^+S,,{\-y^ 


-^) 


y 




i+y'+^ 


- COS 9 


V 



The integrals above are finite and can be easily evalu- 
ated in large momenta limit, p^ 3> 4, introducing the 
dimensionless scale a ~ Am. 

As an example we can evaluate the vacuum polariza- 
tion operator for gi wavelet. For gi wavelet the regular- 
izing function 



where p is dimensionless, i.e. is expressed in units of to,. 



Introducing the notation /3{y) = 



V 
,2?VPii 



^— and using the 



Ay^Sf,^ + By^^^, we get 



substitution yp^y^ 

Yl(A) g2 roo 

__|__ = __(^2p2)y^ dyyFA{mp,mpy) 

. ,/,.((2A-l)y2+i 
X / d9 sm^ 9— 



Pii.pi' 



i2By' - h) 







;32(y)-cos2 



l2„2\ 



Fa{p, q) == cxp ( - A2p2 _ 4^2^- 

Hence for large p^ ^ 4 the integral (36) can be evaluated 
by substitution y^ — t [20]: 



AA) 



,2„2 



"67r2™^ I SaV 



{'■ 



(4aV - a^P^ - 1) + 



8a6p6 



(l_4aV + 2aV) - ^Eii(aV) +Eii(2aV)}- 
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Similarly, the longtitudinal term X^ evaluated with gi 
wavelet in the hmit p^ ^ 4 is equal to 



167r2 



a^p^ 



In the limit of small scales ap <C 1 the equation (37) 
does not depend on p: X^ ex -^. Therefore the whole 
equation (36) is similar to the result obtained by Pauli- 
Villars regularization of the vacuum polarization 



nff^ip) = cM^S^, + [pH^^ - p^p,,)F ( ^ 



m 
V 4m2 ' M 



where M ~> oo is a regularizing mass [42]. The gauge 
invariance is restored if the multiscale diagram (35) is 
integrated over all scales. In this limit the theory can be 
subjected to dimensional regularization [1]. 

c. Vertex function. The one-loop contribution to 
the QED vertex function for the theory with scale- 
dependent matter fields is shown in Fig. 8 below. The 




a^_ s X P 

p,a \k,a 

FIG. 8: One-loop vertex function in scale-dependent QED 

equation, which corresponds to the vertex diagram 
Fig. (8) can be casted in the form 



-le 






{-^e)' 



dH 



it can be seen that the right-hand side of the equation 
(38) can be represented as a linear combination of three 

finite integrals [J^^^ , J^ , J^J) presented in Appendix C, 
in analog to their divergent counterparts in Minkowski 
space [43]. After some algebra the vertex (38) turns to 
be 






g{-pa')g{-qr)g{ka) 



+ 



7(2)1 
(39) 



C. Ward-Takahashi identities 

The Ward-Takahashi identity in spinor electrodynam- 
ics relates the vertex function to the difference of fermion 
propagators: 



qiXtj.{p,q,p + (i) = G '^{p + q)-G ^(p). 



(40) 



where G{p) is the complete fermion propagator. The 
identity (40) is a helpful constraint which ensures the 
gauge invariance of the renormalized QED in any order of 
perturbation theory [44, 45]. The constraint (40) makes 
the perturbation expansion gauge invariant at the pres- 
ence of the gauge fixing terms in the QED generating 
functional. 

In wavelet-based theory, where the fields explicitly de- 
pend on scale the divergence do not appear in the Feyn- 
man diagrams, but the evaluation of integrals in internal 
lines with the integration scales constrained by the mini- 
mal scale of external lines may spoil the gauge invariance 
of the complete propagator. To prevent this the Ward- 
Takahashi identities are required. 

At the absence of gauge- fixing terms in the Lagrangian 
(27), the generating functional 



-gi-pargi-qrygika) ^ ^"^W (2^)^ '"^^" '" ""' 
xG(k ~ fYipD^pFA^P - f)FA(k - f)FA{f). 

The explicit substitution with photon propagator taken 
in Feynman gauge gives 



7„G(p - /)7^ X e-^I^^'^'X] ^ / VAVi;ViPe- ^ d^^iLi4,.4,.A)+JA+^xi:+^i:x) ^ 



le 



p(^) 



g{-pa')g{-qr)g{ka) 



i->ef • '^ 



■la 



-f-m 



(27r)4 '" {p - ff + m2 

— / — Tfl 1 

X7m „: L , ^, 1o.^Fa{p- I)FA{k- DFaU) 



(41) 

with L{ip, if), A) given by the equation (27), would be in- 
variant under the gauge transformations (26) if no source 
term — J d^x{J^A^ + txip + I'lpx) is present. 

In the framework of scale-dependent functions the 
gauge field A^{x) is expressed in terms of its wavelet 
coefficients A^a{b)- 



■ (fc - /)2 + m2 '" /2 



(38) 






b\ , ,,.dad% 

Af^aib) 



Representing the numerator of the latter equation in the 
form 

^/.= la{f-f)lt.{$-f)la- 

m [(|^ - /)7^ + 7^(1^ - /)] + 2mV: 



(with the angular part of wavelet transform (5) dropped 
for simplicity) . In view of linearity of wavelet transform 
we may infer the gauge transform of the scale components 
to have the form 



Ka{x) = A^a{x) + 



dK{x) 
dx„ 
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where 



Kix) 



jjd a' 



]^-g{lL_l]Aiy)d^y 



variation the integrand in the functional integral (41), 
after integration by parts, acquires a multiplicative factor 
e''^, with 



is the scale component of the gauge function (26). That 
is the gauge transform of the abelian gauge field A^a (x) 
is a projection of the (no-scale) gauge field A^(x) onto 
the space of resolution a. 

Since the free Lagrangian of QED is gauge-invariant 
by construction, the derivative of the Ward-Takahashi 
identities turns into evaluation of the functional overage 
of the variation of source and gauge fixing terms under 
inifinitesimal gauge transform 



Sa= I d^'x --d^{d^A^) + d^J^ + eii;x~Xi^) A(a;). 

(42) 

Considering Sa as small we can e '^ « 1 + (Ja and proceed 
with the derivation procedure from 



5 - (Sa) - 



(43) 



SAf^ = d^K, Sip = —leAip, Sip = leKip, 
where A — A(a;) is considered to be small. Under this 



The standard procedure of the variation of action with 
a gauge fixing term [46] with respect to Aa(x) (43) leads 
to the equations [47]: 



J 



q^^^aia3a^{ji),q,P + q)= / ^G^iq^ (P + g)Ma2a3a4 (P + 9, <?, P) - / ^^aiagaa (P + 9, <?, P)^^,^^ (p), (44) 



where 



02 



da2 

02 



^^010203(^1,^2,^3) = (27r) S (fci - fc2 - h)giaiki)g{a2k2)g{a3k3). 
I 



The equation (44) is exactly the wavelet transform of the 
standard Ward-Takahashi identity (40). 



D. QCD example 

Same as in QED we can evaluate the gluon vacuum 
polarization operator using 171 as the basic wavelet. The 
corresponding one- loop diagram is shown in Diagram 45. 




^■'-<LM= (45) 



2^ ^ J (27r)4 P{1+pY 



is the tensor structure of the vacuum polarization dia- 
gram (45) in M^ Euclidean space. A is the minimal scale 
of two external lines. The regularizing function, if calcu- 
lated with gx wavelet, has the form (32): 

Fa{1+P, I) = exp(-2yl2(/ +p)2 _ 2A^P). 

Symmetrizing the loop momenta in equation (45) by 
substitution I = g — | , we obtain 



^^^^.M= -^/■'^^r^^ 



AB,tj.v 



d^q 

(2^ 



FA(p,q)x 



Wq^Qt, - ^Pf,p^ + (5^^(|p2 + 2g2) 



„2 _ Hi 

q 4 



(46) 



where 



N^v{l,p) = I0ltj.liy + 5{lfj,Pu + luPp.) - 2pfj,p^ 

+ {p- lfS^„ + {2p + ifS^, 



For gi wavelet the regularizing function Fj^(p, q) is given 
by equation (32). 

The integral (46) can be easily evaluated in infrared 
limit where ordinary QCD is divergent: 



ni1^:';l(p^o) = -,V^^.f^^^ 



d'^q e 



-AA^q' 



(27r) 



-[5q^g,. + q'^Sfj,ty] 



(47) 



Making use of the isotropy 



d q ^ 2^ q dq, q^^q^ -> Sf_,^ — 



we get 



2 fACD fBDCi 



Qg.f''^''.f 



32 



qdqe 



gg2 fACD fBDCs^^ 

256A^ 
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Similar contribution conies from the ghost loop. 
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V. CONCLUSION 

In this paper we developed a regularization method 
for quantum field theory based on a continuous wavelet 
transform. Regardless significant amount of works de- 
voted to wavelet-based regularization in different quan- 
tum field theory models [11, 12, 48], all those are ba- 
sically the lattice theories. The novelty of the present 
approach, developed by the authors [1, 13, 49], consists 
in using continuous wavelet transform to substitute the 
local fields (j){x) by the scale-dependent fields 4>a{x), de- 
fined as wavelet-coefficients of the physical field. Sub- 
stitution of such fields into the action, supplied by ap- 
propriate causality assumptions and operator ordering 
[1, 2, 28], results in effective regularization of Feynman 
graphs, which makes each internal line decay as an effec- 
tive factor ex e^P , where A is the minimal scale of all 
internal lines, and p is momentum. 

Regulazation factors, that are technically similar to 
our approach, were already known in QCD. They are re- 
lated to the modification of the gluon vacuum state to 
the instanton vacuum, with the parameter A understood 
as the size of the instanton [50, 51]. The difference be- 
tween the instanton vacuum model and our model is that 
the scattered quark fields are local fields in the instanton 
model and only the interaction with instanton vacuum 
is smeared. In our approach the incident particles are 
nonlocal wave packets and only the integration over all 
scales makes the theory local. 

The physics of using scale-dependent fields (j)a{x) in- 
stead of local fields (j){x) lies in the fact, that no physical 
quantity can be measured in a point, but in a region of 
nonzero size a > 0. Thus only the finite resolution pro- 
jections (j)aix) of a quantum field (f) are physically mean- 
ingful. The n-point Green functions for such fields con- 
structed by our method are finite by construction and do 
not require regularization. The gauge invariance of the 
theory results in appropriate Ward-Takahashi identities, 
which are the projections of ordinary Ward-Takahashi 
identities onto finite resolution spaces. 

The practical applications of our approach can be 
found in such physical settings where the separation of 
the field to the components of different scales is physi- 
cally meaningful. Such models have been presented in 
QED calculations of the dependence of the Casimir force 
on the size of displacement in measurement [49] , and also 
in application of quantum field theory methods to the 
calculation of correlations of the turbulent velocity fiuc- 
tuations of different scales [52]. We strongly hope that, 
regardless the yet unsolved problem of deriving renor- 
malization group equation in continuous limit of wavelet- 
based theory, this method can be also applied for QCD 
calculations, where it was originally proposed [11]. 
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Appendix A: Dirac 7-niatrices in Euclidean space 



7m7p == -4, 7p|*7p = 2|* 



(Al) 
(A2) 



Slashed vectors denote convolution with Dirac gamma- 
matrices ^ = 7pfcp, ^^ = — fc^ 

Appendix B: Feynman rules in Euclidean space 

Photon propagator is taken in Feynman gauge: 



D{k) = :g 



Fermion propagator 

Gf{p) = 



-I f - m 

— i- 



j + m p^ + m? 
Electron-fermion vertex: 

Besides that each fermion vertex results in extra sign 
of the whole diagram. 

Appendix C: Functions and Integrals 

Exponential integral of the first type 



Eii(z) = / dx 



Integrals for angle integration in Euclidean Green func- 
tions [20]: 

r ,^ sin^ e cos^''^ e 
io{y) = 7r(i - v/i-r'(y)), (Ci) 
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The constants A, B for the vacuum polarization diagram 
(35) are given by 4A + B = 1, A + B = /i//o, from 
where we get 



A 



1 



1 



+ I-\y)- l3\y)^B 



\-^I,^i,)40^iy)■ 



3 6 " '•" 3 
Integrals in one-loop fermion-photon vertex 

d^f FA{p-f)FA{k-f)FA{f) 



J(0) 



j(l' 

"^M 



7(2) _ 



(27r)4 [{p - /)2 + to2][(/c - /)2 + m2]/2 

c^V UFa{p - .f)FA{k - f)FA{.f) 
(27r)4 [(p - /)2 + m^][{k - /)2 + m2]/2 

d^I U.UFA{p-f)FA{k-f)FA{f) 
(27r)4 [(p - /)2 + to2][(A: - /)2 + m2]/2 



(C2) 
(C3) 
(C4) 
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